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PREFACE 


Tbe  work  described  in  this  Memoraivlum  was  pursued  in  connection 
with  a  study  of  radar  reflection  and  absorption  characteristics  of 
the  wake  of  a  re-entry  vehicle.  The  study  is  one  of  a  series  related 
to  the  ICBM  discrimination  problem  currently  being  pursued  at  RAND 
under  the  sponsorship  of  the  Advanced  Research  Projects  Agency. 


SUMMARY 


The  problem  of  the  absorption  and  reflection  of  electromagnetic 
energy  from  the  wake  of  a  body  moving  at  hypersonic  speed  througli 
the  atmosphere  has  engaged  the  attention  of  a  nvunber  of  researchers 
in  the  defense  community.  In  connection  with  that  problem,  this 
paper  presents  a  simplified  analytical  model  for  investigating  such 
wake  characteristics  as  electron  density,  collision  frequency,  velocity, 
and  temperature.  Thus,  the  model  fxirnishes  the  information  which  is 
necessary  for  studying  the  electromagnetic  scattering  properties  of 
the  wake. 

The  model  is  based  op  the  work  of  Paul  lykoudis,  and  assumes 
inviscid  laminar  flow  about  a  blunt  body.  A  simplified  expression 
for  the  enthalpy  in  the  wake  is  derived  from  lykoudis'  universal 
solution.  The  enthalpy  is  then  translated  into  the  desired  wake 
characteristics  by  using  the  equilibrium  properties  of  air.  For 
simplicity  the  air  is  assumed  to  be  in  local  equilibrium.  This 
represents  an  approximation  since  the  air  is  strictly  in  local 
equilibrium  only  at  lower  altitudes  where  the  wake  is  turbulent. 

A  particular  set  of  re-entry  conditions  is  considered  as  eui 
example,  and  an  analytical  e:^resslon  for  the  electron  density  is 
derived.  This  expression  is  compared  with  the  electron  density 
contours  derived  in  previous  studies  elsewhere,  and  the  differences 
are  discussed.  Finally  the  general  structure  of  the  wake  is  illustrated 
by  plotting  electron  density  contours. 
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I.  nmoDucTioN 


Recently  there  has  been  a  great  deal  of  Interest  In  the  problem 
of  absorption  and  reflection  of  electromagnetic  energy  from  the  ionized 
trail  or  wake  left  by  a  body  moving  at  hypersonic  velocities.  In 
order  to  determine  the  electromagnetic  properties  of  the  wake,  it  is 
necessary  to  know  in  some  detail  the  structure  of  the  we^e  in  terms 
of  the  In^rtant  electromagnetic  peo-ameters  such  as  electron  density, 
collision  frequency,  velocity,  temperature,  etc. 

The  structure  of  the  wake  has  been  investigated  by  many  authors. 

(1  2) 

Using  numerical  techniques,  Feldman'  *  '  was  one  of  the  first  to 
examine  the  wake  structure.  Recently,  tykoudls^^^  has  investigated 
the  problem  using  an  analytic  approach  throu^out.  Any  of  the  above 
references  will  provide  the  reader  with  an  excellent  introduction  to 
the  fundamental  nature  of  the  wake.  It  is  the  piuT>ose  of  this  note 
to  transform  the  analytic  results  of  lykoudls,  which  are  in  terms 
of  the  enthalpy,  into  the  above  electromagnetic  parameters.  At  the 
same  time  an  attempt  will  be  made  to  collect  and  reorganize  the 
necessary  collateral  data. 

In  Section  II  a  brief  general  discussion  of  the  wake  will  be 
given.  Section  II  will  Introduce  the  following  sections,  and  will 
discuss  the  basic  assumptions  and  limitations  of  the  model.  In 
Section  III,  lykoudls'  results  for  the  conduction-controlled  part  of  the 
wake  are  presented,  and  in  Section  IV  we  provide  the  necessary  data  for 
translating  the  lykoudls  model  into  the  important  electromagnetic 
pareuneters,  and  an  approximate  analytic  model  for  the  electrcnagnetlc 
properties  is  derived. 
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II.  GENERAL  DISCUSSION  OF  THE  WAKE 

TSie  structure  of  the  hypersonic  wake  behind  a  blunt  body  is 
shown  in  Fig.  1.  For  a  detailed  discussion  of  the  structure  the  reader 
is  referred  to  Ref.  3.  Referring  to  Fig.  1,  we  may  briefly  describe  the 
flow  pattern  as  follows.  The  free-stream  air  is  irreversibly  compressed 
while  passing  throvigh  the  bow  shock.  The  air  then  expands  and  is  re¬ 
compressed  in  the  trailing  shock.  This  is  followed  by  further  e^qansion 
and  at  some  station  (x  a  x^)  the  air  has  expeuided  to  the  free  stream 
pressure  p^. 

The  portion  of  the  flow  which  enters  the  boundary  layer  forms 
the  viscous  wake.  Beyond  the  station  x  =  x^  the  wake  is  characterized 
by  the  diffusion  of  energy  out  of  the  central  viscous  part  of  the  wake. 
In  a  Laminar  wake  this  transport  takes  place  by  molecular  diffusion; 
in  a  turbulent  wake  the  process  takes  place  by  txarbulent  diffusion  and 
mixing. 

Our  dlscvission  will  be  restricted  to  blunt  bodies.  For  slender 
bodies  the  shock  is  weaker  and  the  viscous  wake  contains  most  of  the 
momentum  loss  and  hl^-teraperature  air.  As  a  resiilt  the  electromag¬ 
netic  parauneters  depend  critically  on  the  boundary  layer  at  the 
trailing  edge  of  the  body,  and  hence  on  the  afterbody  shape.  For 
blunt  bodies  it  has  been  shown  that  the  exact  shape  of  the  afterbody 
does  not  appreciably  affect  the  nature  of  the  wake.  Ibis  occiurs 
beca\i8e  the  Invlseld  flow,  which  is  relatively  vinaffected  by  the  after- 

I 

body  and  trailing  shock,  contains  most  of  the  momentimi  defect  by 
means  of  the  strong  bow  shock,  and  hence  the  bow  shock  controls  the 
behavior  of  the  wake  in  this  case.  We  shall  asstune  that  this  is  the 


case  emd  thus  we  shall  Ignore  the  contribution  of  the  trailing  shock 
to  the  pressure  and  the  velocity.  Neglecting  the  afterbody  and  neck 
shock  is  equivalent  to  assvuning  that  the  wake  may  be  computed  from 
the  nature  of  the  bow  shock  alone. 

We  shall  also  neglect  the  effects  of  the  boundary  layer  contri¬ 
butions  to  the  viscous  wake  suid  consider  instead  that  the  far  laminar 
w8dce  is  a  gas  cylinder  with  a  given  initial  temperatiure  distribution 
which  cools  by  conduction  suid  convection.  IRils  assumption  is  tantamount 
to  the  neglect  of  turbulence,  where  the  principal  radial  cooling 
mechanism  is  turbvilent  diffusion.  It  appears  that  this  assumption  is 
reasonable  for  a  significant  range  of  flight  conditions,  in  peortlcular 
for  blunt  bodies  at  altitudes  above  180,000  ft  and  below  250,000  ft. 

At  altitudes  above  250,000  ft  we  violate  our  assumption  that  the  flow 
is  Invlscid;  emd  at  altitudes  below  about  180,000  ft  the  flow  is  likely 
to  be  turbulent,  and  thus  we  violate  our  assiuoptlon  that  the  flow  is 
laminar.  A  detailed  discussion  of  the  structure  of  the  turbulent  wake 
will  be  found  in  Ref.  U. 

The  effects  of  turbulence  on  the  electromagnetic  properties  of 
the  wake  are  two-fold.  First,  tiirbulence  tends  to  cool  the  wake  much 
more  rapidly,  and  second,  turbxilence  Introduces  local  fluctuations  in 
the  parameters.  This  complicates  the  scattering  problem  by  aerody- 
namlcally  "roiighenlng"  the  surface  of  the  wake.  This  causes  inter¬ 
ference  of  the  incoming  and  reflected  electromagnetic  waves,  and  in 
general,  significantly  changes  the  (unplitude  of  the  scattered  energy. 

The  meignltude  of  the  tvirbulent  effect  is,  of  course,  a  function  of  the 
scale  of  the  turbulence  as  well  as  the  rcuiar  frequency  and  mean 
electron  density. 
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The  normal  procedure  for  analyzing  a  wake  is  to  divide  the  wake 
into  two  distinct  regions,  one  of  which  is  controlled  by  expansion  of 
the  invlscld  flow,  and  the  other  determined  by  an  energy  balance  of 
axial  convection  and  radial  conduction  from  the  inner  wake.  Hie  ex¬ 
pansion  region  starts  (when  the  neck  shock  is  neglected)  at  the  bow 
wave  and  continues  until  the  diffusive  or  molecular  processes  become 
important,  somewhere  in  the  vicinity  where  the  pressure  approaches 
ambient  values.  For  convenience,  we  will  assume  that  the  expansion 
region  Includes  all  the  flow  from  the  bow  shock  up  the  point  x  =  x^ 
where  p  =  p  •  In  this  region  the  flow  is  asstuned  to  be  Invlscld  and 
adiabatic . 

The  second  region  is  called  the  conduction  region  and  is  the 
remainder  of  the  wake.  In  this  region  the  only  important  mechanism  is 
assumed  to  be  the  cooling  of  the  wake  by  thermal  conduction. 

In  some  cases  a  significant  amount  of  conduction  cooling  takes 
place  before  the  pressure  reaches  ambient  values,  at  station  x^.  For 
our  two-region  model  to  be  valid,  this  conduction  cooling  must  be 
small  coii5«red  to  the  exapnslon  process.  A  detailed  study  of  the 

(5) 

relative  magnitudes  of  these  two  processes  has  been  made  by  lykoudis.'" 
His  results  indicate  that  for  body  radii  greater  than  where 


r 


min 


(1) 


the  two-region  model  is  adequate. 

For  typical  re-entry  conditions  the  Reynolds  number  is  high  and 
Eq.  (1)  indicates  that  our  two-region  model  would  be  adeqviate  for 
bodies  larger  than  1  cm  in  radius. 
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For  altitudes  above  200,000  ft  the  flow  may  be  conbi„.  red  "frozen" 
for  bodies  which  have  a  nose  radius  of  1  ft  or  greater,  and  the  composi¬ 
tion  of  the  wake  air  may  be  considered  to  be  the  same  as  the  composition 
of  equilibrium  air  under  stagnation  conditions.  On  the  other  hand,  for 
altitudes  below  100,000  ft  the  relaxation  rates  are  sufficiently  fast 
so  that  the  air  may  be  considered  to  be  in  thermodynamic  equilibrium 
with  the  local  conditions.  In  order  to  simplify  the  calculations  we 
shall  assume  that  the  air  is  in  thermodynamic  equilibrium  with  the 
local  conditions  at  all  times. 

The  effect  of  ablation  products  on  the  electron  density  and  other 
electromagnetic  parameters  can  be  very  important.  With  some  reluctance 
we  will  neglect  the  ablation  products  in  ovir  analysis.  There  are  two 
reasons  for  this.  First,  very  little  information  is  available  on  the 
equilibrium  characteristics  of  air  containing  typical  concentrations 
of  ablation  products  and  second,  since  we  have  neglected  the  boundary 
layer  there  is  no  appropriate  way  to  introduce  the  ablation  products 
into  the  problem. 

In  Ref.  3,  it  is  pointed  out  that  the  length  of  the  expf’sion 
region  is  small  compared  to  the  length  of  the  conduction  region.  The 
length  of  the  expansion  region  may  be  estimated  from 


X 

P 


r 

o 

ITT 


(2) 


Since  the  expansion  region  is  small  we  will  consider  only  the  conduction 
part  of  the  wake.  If  the  reader  wishes  to  Include  the  expansion  region 
in  his  model  he  may  use  the  results  presented  in  Ref.  3> 

In  the  next  section  we  will  present  lykoudls*  results  for  the 


conduction  region. 
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III.  SIMPLIFIED  MODEL  OF  TtCE  CONDUCTION  REGION  OF  THE  WAKE 

In  this  section  we  present  ly’toudis'  universal  solution  for  the 
enthalpy  in  the  conduction  region.  The  first  part  of  this  solution 
determines  the  enthalpy  profile  at  the  station  x  =  x^.  lykoudis 
obtained  the  following  result  from  the  conservation  of  total  enthalpy 


h(x  .r) 


Vr; 


(5) 


where  we  have  introduced  the  Howarth  radius  R  which  is  related  to  the 
actual  radius  r  by 


2  I*  ^ 

The  term  is  an  effective  ratio  of  specific  heats  which  is  approximately 
1.2  for  hypervelocity  flight  conditions.  The  parameter  a  and  the 
drag  coefficient  determine  the  shape  of  the  bow  shock.  Reference 

4  gives  the  following  relation  for  the  bow  shock 


R  = 


(a  C_) 

i) 


1/4 


a  is  approximately  2.6  and  f  is  defined  below; 


f 


(5) 


(6) 
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f  varies  betveen  0.77  and  0.85  as  M  varies  Isetween  15  and  55. 

2) 

Feldman'  ’  '  has  numerically  computed  the  enthalpy  profile  at 

(5) 

X  =  Xp  by  using  the  method  of  characteristics.  lykoudis'"^'  malces  a 
comparison  between  Feldman's  results  and  Eq.  (5).  The  two  results 
agree  very  well.  lykoudis^^^  and  Feldman^ have  also  compared  a 
gaussian  profile  (to  be  considered  later)  with  Feldman's  results. 

The  gavisaian  profile  agrees  with  Feldman's  resvilts  for  small  R,  but 
diverges  for  R  >  r^,  where  Eq.  (5)  becomes  a  better  approximation. 

V/e  note  that  Eq.  (5)  may  be  approximated  by 


h(Xp,r) 

h(Xp,o) 


1  + 


C  r^ 
D  o 


(7) 


A  solution  for  the  enthalpy  in  the  conduction  region  may  be  found 

(assuming  constant  thermal  conductivity)  for  the  boundary  condition 

given  in  Eq.  (5).  This  solution  is  given  in  terms  of  a  complicated 

Fourler-Bessel  IntegreQ,.  In  order  to  simplify  the  problem,  let  us 

assume,  following  Feldman^^’^^  and  lykoudls,^^^  that  the  initial 

distribution  at  x  =  x^  is  gaussian. 

P 


(8) 


Neither  Eq.  (3)  nor  Eq.  (6)  are  valid  for  RhRg,  where  Rg  Is  the 


9 


radius  of  the  bow  shock  at  x  =  x  .  For  R  ^  R_  the  enthalpy  Is,  of 

P  ° 

course,  equal  to  the  ambient  enthalpy,  h^. 

A  comparison  of  Eq.  (8)  with  Eq.  (7)  Indicates  that  the  two  dis¬ 
tributions  match  very  well  for  small  R,  and  as  long  as  we  restrict  our 
solutions  to  small  R  no  serious  errors  will  arise. 

The  use  of  Fq.  (8)  has  the  advantage  that  the  energy  equation  in 
the  conduction  region  can  be  solved  In  closed  form  even  when  the  thermal 
conductivity  Is  variable.  It  can  be  shown  from  Ref.  (5)  that  the 
enthalpy  distribution  for  regions  downstream  of  x  »  x^  Is  given  by 


(1  +  X) 


'0.8 


(1  +  X) 


o 


S  ^o 


(9) 


where 

1/U 

I 

^  "  P-  C  r^ 

0.554?  —  0 


In  cgs  vuilts. 

Now  for  small  R  we  can  replace  the  exponential  In  Eq.  (9)  toy: 


(1  +  Cp  rl 


1  + 


(10) 


(1  +  X)--  Cj, 
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Therefore  Eii.  (9)  may  be  written  as 


/  nO-8 
(1  +  X)  + 


C  r^ 
D  o 


(11) 


We  note  In  passing  that  the  substitution  we  have  Just  made  Is  not 
likely  to  Introduce  additional  errors  Into  our  computation,  but  to 
the  contrary  It  is  likely  to  improve  our  answer  since  Eq.  (ll)  satisfies 
the  more  accurate  boundary  condition  Eq.  (7). 

We  summarize  our  universal  solution  as  follows; 


h(x,r) 


R  T 


0 


h(Xp,o) 


R  T 

o 

• - - - - - 

(1  +  x)  +  ■  '  A 

Cn  ^ 

D  o 

(12) 


/h(Xp,o)Y/'* 

V  R  ^0  /  ^ 

0.5545  —  U  ^0^ 

Po  • 


(13) 


2 

r 


(14) 


In  the  next  section  we  will  discuss  the  numerical  computation  of 
h  and  its  translation  into  the  electromagnetic  parameters. 
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IV.  DETmMINATION  OF  THE  WAKE  PROProTIES  FRCTl  THE  SIMPLIFIED  MODEL 

In  this  section  we  v;ill  provide  the  data  for  determining  the 
parameter  profiles  for  a  given  set  of  re-entry  conditions.  We  assume 
that  the  vehicle  velocity  (free  stream  velocity  U^),  the  altitude  H, 
the  drag  coefficient  and  the  body  radius  r^  are  given.  The  terms 
and  H  are  taken  from  the  trajectory  of  the  body;  examples  are  shown 
In  Fig.  2. 

The  first  step  Is  to  translate  our  altitude  data  into  the 
required  ambient  or  free  stream  thermodynamic  parsuneters.  To  do  this 
we  will  assume  the  ARDC  model  atmosphere. The  thermodynamic  vari¬ 
ables  are  plotted  as  a  function  of  altitude  In  Figs.  7  through  12  In 
Appendix  A.  For  simplicity  we  will  assume  the  ambient  conditions  do 
not  change  appreciably  over  that  portion  of  the  wake  in  which  we  are 
Interested,  and  thus  may  be  taken  as  constant.  The  range  of  ambient 
conditions  should  be  checked  and  an  appropriate  average  value  chosen. 
If  a  more  acciirate  solution  is  desired,  the  wal^e  may  be  divided  up 
into  a  number  of  sections  where  the  ambient  conditions  of  each  section 
may  be  considered  constant. 

Our  next  step  will  be  to  compute  h(x,r)/RT^  from  T3qs.  (12), 

(15),  and  (14).  The  only  quantity  appearing  in  these  expressions 

which  cannot  be  read  directly  from  Figs.  ^  through  12  is  h(x  ,o)/RT  . 

(7)  P  o 

This  quantity  may  be  computed  from  shock  tables'  '  or  from  the 
approximate  relation  below,  which  is  derived  in  Ref.  4  from  the  energy 
equation  assuming  that  f  =0.3. 
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h(x^,o) 
R  T 


a  Jk.  (y  -1) 


R  T 


(15) 


'fhere  M  =  —  .  The  term  h(x,r)/RT  is  now  completely  determined. 

We  next  turn  to  the  problem  of  determining  the  thermodynamic  and 

h  ^o»  , 

electromagnetic  parameters  from  —  for  a  given  pressure  —  (i.e., 

A  T  p 

O 

altitude).  First  let  us  consider  the  average  velocltjr^of  the  air  in 
the  wake.  For  a  constant  stagnation  enthalpy  in  the  radial  direction 
(Prandtl  number  unity)  the  conservation  of  energy  implies  that 


1 


(16) 


For  h/h^  »  1  and  U/u^  «  1  the  above  expression  reduces  to 


,  U(x,r) 

u. 


2  .  Ik 


which  can  be  compared  with  its  value  at  x  =  x  ,  r  =  0  to  yield 


_  u(.x^r2 ,  r  1 . 

L  U,  -J 


(17) 
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We  recall  that 


U(x^,o) 

~U 


f  w  0.8 


(18) 


so  therefore  the  air  velocity  with  respect  to  a  stationary  atmosphere 
is  given  hy 


V  .  0.  -  U(x,r)  .  0.2  U.  (191 

P 

The  remaining  parameters  are  presented  graphically  as  a  function 

of  h/RT  for  various  values  of  log  —  under  the  assvunption  that  the 

^o 

air  is  locally  in  equilhrium.  The  S,  p,  T  and  z  curves  were  obtained 

from  Blackwell, the  N  curves  from  Gilmore, and  the  u  curves 
'  e  e 

from  Romig.  These  curves  are  shown  in  Figs.  15  throu^  19>  in 
Appendix  B. 

In  some  applications  one  prefers  to  have  the  structure  of  the 
wake  presented  in  graphical  form;  for  example  electron  density  con- 
tovirs.  In  this  case,  the  most  efficient  way  to  obtain  the  graphical 
representation  is  simply  to  plot  it  point -by- point  from  the  graphs 
presented  in  Appendix  B.  However,  in  most  applications  the  structure 
of  the  wake  serves  as  a  starting  point  for  the  analysis  of  a  particular 
problem;  for  example,  the  computation  of  the  radar  cross  section. 

In  this  case  it  is  much  more  convenient  to  have  the  structure  of  the 
wake  presented  in  analytic  form.  In  this  section  we  will  illustrate 
the  technique  of  determining  an  approximate  analytic  representation 
by  considering  a  specific  example. 

‘N’ 

Private  communication. 
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The  input  conditions  for  the  specific  case  to  be  considered  are 
given  below; 


U  =  25,000  ft/sec  =  7*62  X  10^  cm/sec 


H  =  100,000  ft 


(20) 

(21) 


The  terms  Cj^  and  r^  will  be  retained  as  parameters  and  r^  is  to  be 
given  In  centimeters.  This  altitude  was  selected  so  that  the  results 
could  be  compared  with  Feldman,  however  the  wake  is  likely  to  be 
turbulent  at  this  altitude. 

Using  these  input  conditions  we  readily  obtain  the  following 
results: 


From  Fig.  7 


CB 

R  T, 


=  2.78 


(22) 


From  Fig.  8 


P, 


1.28  X  10 


-2 


(25) 


From  Fig.  10 


log  —  =  -  1.96 
•^o 


i2k) 


From  Fig.  11 


a  =  99lf  ft/sec 


(25) 


Therefore 


M  =  25.2 


(26) 
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From  Eq.  (15) 


h(x  .o) 


From  shock  tables 


h(x^,o) 


V7e  shall  use  the  value  obtained  from  the  shock  tables,  since  Eq.  (15) 
gives  only  an  approximate  value  for  h.  For  this  partictilar  case 
Eq.  (15)  gives  a  value  which  is  about  15  per  cent  below  the  correct 
value. 


From  Eq.  (I3) 


X  =  6.64  X  10 


.4 


D  o 


Prom  Eq.  (2) 


140  r 


From  Eq.  (12) 


h(x,r) 
R  T 

o 


(1  .  x)°-®  .  ^ 


r 

D  o 


Since  we  will  only  be  interested  in  the  structure  of  the  wake 

vhere  "  2  JO,  or  for  electron  densities  greater  than  10^ 

R  T 

electrons /cm' ,  the  overall  range  of  ^  to  be  considered  is 

■^o 


50  s 


s  150 


(52) 
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Next  we  will  determine  the  density  as  a  function  of  -  _  . 

R  Tq 

This  relation  is  given  graphically  in  Fig.  14.  The  density  ciirve  can 
be  approximated  by  a  straight  line  to  within  a  factor  of  1.1  for  the 
range  of  values  given  in  E^.  (52).  The  straight  line  on  log-log  paper 
gives  the  following  power  law  relation: 


or 


From  Eq.  (l4) 


f  ■  <irV) 


•0.6 


(33) 


£_ 

P, 


0.78  (^)‘°'^ 

00  O 


vf 


2R  dR 


0.78 


142 


1  -0.6 


(1 .  x)°-®  *  -£ 


Cn 

D  o 


(35) 


Integrating  Eq.  (35)  gives 


,,  „^0.8  R 

(1  +  X)  +  - 5- 

D  0 


(1  +  x)°’^^  +  0.016  — ^ 

CT^  r 
D  o 


2.5 


(36) 
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Substituting  Eqs.  (29)  and  {36)  into  Eq.  (3I)  gives 


R  T 

o 


IU2 


1  +  6.64  X  10' 


0:32 


I  2.5 


+  0.016 


D  0 


(37) 


The  other  quantities  may  be  determined  as  a  fi  notion  of 

o 

by  empirically  fitting  the  appropriate  curve  for  the  given  range 

of  ^  ■  .  Since  the  other  quantities  depend  only  on  the  enthalpy, 

o 

the  constant  enthalpy  contours  may  be  interpreted  as  contours  of  the 
other  quantities  also. 

Let  us  consider  the  electron  density  relation  as  an  illustrative 

example  of  the  technique.  The  electron  density  relation  is  shown  in 

Fig.  18.  The  curve  may  be  fitted  to  within  a  factor  of  5  by  a  stral^t 

line  for  the  given  range  of  ^  . 

■‘'o 


N 

e 


2  X  10'®  (^)^'^ 


electrons 


cm 


Substituting  Eq.  (37)  into  (38)  gives 


(38) 


2  X  10 


12 


e  r 


1  +  6.64  X  10 


-4 


+  0.016 


D  o 


electrons 


(39) 


cm' 
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The  other  quantities  such  as  temperature,  collision  frequency,  etc., 
may  be  found  in  a  similar  fashion. 

For  other  re-entry  conditions  one  obtains  expressions  for  the 
electron  density  which  are  very  similar  in  form  to  Eq.  (59).  For 
example,  H  =  100,000  ft  and  U  =  20,000  ft/sec  gives 


N 


2.k  X  10 


10 


( 


1  +  7.21  X  10 


-4 


[x  - 


C  r^ 
D  o 


) 


(5732" 


+  0.021 


C  r^ 
D  o 


■55^2 


(Uo) 


In  Figs.  3  and  4  we  compare  Eqs.  (59)  and  (40)  with  the  numerical 

(2) 

results  of  Feldman.  Figure  3  shows  the  radial  variation  in  the 
electron  density  at  x  »  x  and  Fig.  4  shows  the  axial  variation  in 
the  electron  density  for  r  =  o. 

Figure  5  indicates  that  the  results  are  in  close  agreement 
for  small  r,  but  tend  to  diverge  for  large  r.  This  divergence  may 
be  traced  to  the  assumption  given  in  Eq.  (10).  Feldman  retains  the 
exponentlaL  form  corresponding  to  the  left-hand  side  of  Eq.  (10),  while 
we  prefer  to  use  the  inverse  form  corresponding  to  the  right-hand 
side  of  Eq.  (10).  Since  the  exponential  form  falls  off  more  rapidly 
with  increasing  r  than  the  inverse  form  does,  the  simplified  model 
will  tend  to  give  a  somewhat  broader  wake. 

Figure  4  indicates  that  there  is  a  basic  difference  between 
lykoudls'  and  Feldman's  expressions  for  the  enthalpy  as  a  function  of 
X.  In  terms  of  electron  density  the  two  solutions  differ  by  a  factor 
of  about  10,  however,  this  corresponds  to  only  a  small  10  per  cent 


Electron  density,  N^,  electrons/cm 


20 


Radial  distonce,  - 

vCo  ro 


Fig.  3 — Radial  electron  distribution 
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AkIoI  distance,  — - y  x  10"’ cm'' 

Cd  Tp* 

Fig.  4  — Axial  electron  distribution 
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difference  In  terms  of  enthalpy.  Thus  In  terms  of  enthalpy  the  two 

solutions  are  almost  the  same.  Indeed,  it  appears  that  the  uncertainties 

introduced  by  the  simplifying  assumptions  which  have  been  made  in  both 

analyses  will  overshadow  the  discrepancies  between  the  two  results. 

The  two  results  are  in  very  close  agreement  for  the  "hot"  portion 

of  the  wake  (l.e.,  small  r,  small  x).  The  "hot"  portion  is  often  the 

region  of  most  Interest.  In  the  outside  or  cooler  part  of  the  wake 

the  simplified  model  gives  solutions  for  the  walce  which  are  somewhat 

broader  and  shorter  than  the  solutions  given  by  Feldman. 

Finally,  the  nattire  of  the  wake  is  shown  in  Figs.  5  and  6. 

Figure  5  is  a  plot  of  Eq.  (39)  and  Fig.  6  is  a  plot  of  Eq.  (40). 

The  axial  distance  is,  of  course,  greatly  compressed  in  comparison 

with  the  radial  distance.  For  a  one  meter  body  the  trail  is  about 

4 

5  X  10  longer  than  it  is  wide. 
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Electron  density  profiles  (Cose 
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V.  COMCIUSIONS 

The  simplified  model  discussed  in  this  paper  is  applicable  to 
o'  at  portion  of  the  re-entry  trajectory  where  the  flow  may  be  con¬ 
sidered  to  be  inviscid  and  laminar.  This  corresponds  roughly  to  a 
range  of  altitudes  from  180,000  to  250,000  ft.  The  viscous  effects 
become  Important  above  250,000  ft  and  the  flow  is  likely  to  be 

turbvilent  below  180,000  ft.  In  this  range  of  altitudes  the  model  gives 

(2) 

results  which  agree  with  the  numerical  results  of  Feldman. 

In  computing  the  electromagnetic  characteristics  of  the  wake, 
the  air  was  assumed  to  be  in  equlllbrluin  with  the  local  conditions. 

This  assumption  is  not  valid  for  altitudes  much  higher  than  100,000  ft. 
Assuming  local  equilibrium  at  altitudes  above  100,000  ft  will  provide 
a  lower  limit  for  the  actual  electron  density.  At  altitudes  above 
200,000  ft  we  may  assume  the  flow  to  'be  "frozen,”  emd  thus  cooipute 
the  electromagnetic  characteristics  from  the  stagnation  conditions. 

For  the  intermediate  or  transitional  altitudes  the  two  limiting 
cases  will  serve  as  upper  and  lower  bounds  for  the  solution. 
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Appendix  A 
(Figures  7-12) 


(6) 


ARDC  MODEL  ATMOSPHIEE 


Enthalpy,  hco/RTo 
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Altitude,  H,  Hm 

Fig. 7 —  Enthalpy 
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Fig.8  —  Density 
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Fig.  10— Log  pressures 


Speed  of  sound, Om/a 
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0  10  20  30  40  50  60  70  80  90 
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Fig.  1 1 — Speed  of  sound 
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Altitude, H,  km 

Fig. 12—  Temperature 
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Appendix  B 
(Figures  15-19) 


PROPERTIES  OP  EQUILIBRIUM  AIR 


(8,9) 


^0 « 1.288  X  10"®  gm/cm® 


100 

Enthalpy  h  /  RTq 


Fig.  15—  Log.  density 


Temperoture,  T/Tq 
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Fig.  16  —  Temperature 
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Entholpy,  h/RTo 


Fig.  17  —  Compressibility  factor 
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Fig.  18a  —  Electron  density 


4o 


10  100  1000 


Entholpy,  h/RT{, 


Fig.  18  b  — Electron  density 
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